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We study the effects of weak disorder on a Wigner crystal in a magnetic field. We show that an 
elastic description of the pinned Wigner crystal provides an excellent framework to obtain most of 
the physically relevant observables. Using such a description, we compute the static and dynamical 
properties. We find that, akin to the Bragg glass phase, a good degree of translational order survives 
(up to a large lengthscale in d = 2, infinite in d = 3). Using a gaussian variational method, we obtain 
the full frequency dependence of the conductivity tensor. The zero temperature Hall resistivity is 
independent of frequency and remains unaffected by disorder at its classical value. We show that the 
characteristic features of the conductivity in the pinned Wigner crystal are dramatically different 
from those arising from the naive extrapolations of Fukuyama-Lee type theories for charge density 
waves. We determine the relevant scales and find that the physical properties depend crucially on 
whether the disorder correlation length is larger than the cyclotron length or not. We analyse, 
in particular, the magnetic field and density dependence of the optical conductivity. Within our 
approach the pinning frequency can increase with increasing magnetic field and varies as n~^^^ with 
the density n. We compare our predictions with recent experiments on transport in two dimensional 
electron gases under strong magnetic fields. Our theory allows for a consistent interpretation of these 
experiments in terms of a pinned WC. 
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I. INTRODUCTION 

Amongst the various effects of strong interactions in 
electronic systems, the possibility of cpystallization of 
electrons, initially predicted by WignerEJ, is one of the 
most exciting. For densities lower than a critical den- 
sity ric, the Coulomb potential energy dominates over 
the kinetic energy and one expects the formation of a 
Wigner crystal (WC), where the electrons occupy the 
sites of a lattice (triangular in two dimensions). A mea- 
sure of this critical density is the dimensionless parame- 
ter Ts defined as the ratio of the Coulomb to Fermi en- 
ergies. Unfortunately, crystallization requires extremely 
low densities which are hard to obtain for three dimen- 
sional systems. The situation is much better for a twi 
dimensional electron gas where Monte-Carlo simulations! 
have shown that the formation of a WC requires a den- 
sity corresponding to > 37. First, it is easier to ac- 
cess such densities and second, the value of ric can be 
increased by the application of a strong magnetic field 
applied perpendicular to the plane of the 2DEG. Two 
dimensional electron gases under strong magnetic fields 
are thus the prime candidate systems in which to observe 
the phenomenon of Wigner crystallization. 

Indeed an insulatingpstatei was observed in experiments 
on monolayer systemg3a'B'El. This insulating state was 
characterized by a diverging diagonal resistivity pxx for 
temperature T ^ and activated behavior for finite T . 



Moreover, the Hall resistivity pxy was found to be tem- 
perature independent and nearly equal to its classic 
value. For the pure systems, approximate calculationsB 
have shown that the WC becomes the lowest energy state 
when the filling factor v which is the ratio of the density 
to the magnetic field equals v <l/h for Ga-As electron 
systems and around v — 1/3 for the hole like systems. 
So it is quite reasonable to interpret the insulating state 
as a Wigner crystal pinned by the incipient disorder in 
the sample. This is consistent with the diverging lin- 
ear resistivity and with measurements of nonlinear I — V 
characteristics. Nonetheless, given the absence of direct 
imaging of the system, the unambiguous identification of 
this phase as a WC was the subject of intense debate. 
Also, the theoretical analysis was complicated by vari- 
ous important factors compared to the simple case pro- 
posed by Wigner: (i) in these systems in strong magnetic 
fields, the phase diagram is very rich, since the interplay 
of interaction and disorder can lead to phases such as 
the integer and fractional quantum Hall effect. The non- 
crystalline phase is not a free electron gas but a strongly 
correlated phase. It is therefore, not obvious that one 
does not need a description anchored in the quantuna 
Hall physics to describe the insulating phase as wellLl; 
(ii) disorder modifies drastically the physical properties 
of crystalline phases. This makes the determination of 
the transport properties which is our only probe of the 
WC crystal so far, much more difficult to compute. 
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On the theoretical side, the question of how to describe 
such pinned crystaUine phases is still open. Our under- 
standing of the effect of disorder on such elastic strjLfe. 
tures stemmed from the pioneering works p£ Larkintj'tll 
for vortices and Fukuyama and Lee (FL)ll3 for charge 
density waves, which predict that the crystal gets pinned 
by disorder and perfect translational order is lost. As a 
result of the pinning, the a.c. transport develops a peak 
at a disorder dependent pinning frequency ujp. However, 
due to the qualitative nature of these theories, very little 
was known beyond these qualitative aspects of a pinned 
system. This is a tremendous handicap for the WC since, 
contrarily to CDW, the problem of a WC involves many 
length scales. Moreover, the loss of translational order 
and the possibility of defects led to doubts about the va- 
lidity of an el astir theory to study the transport. Recent 
experimental3'Eflli§Ej have started to obtain the full fre- 
quency dependent conductivity, and the magnetic field 
and density dependence of observables such as the pin- 
ning frequency. The fact that these experiments did not 
comply with simple FL type expectations, cast doubts on 
the interpretation of this experimental phase as a pinned 
WC. This clamors for a quantitative theory of pinned 
WC that goes beyond the FL approximations, and with 
which experiments can be compared to. 

Fortunately, the recent progress made in tljp iijjHpr- 
standing of such disordered elastic structuresllZloO'EHl 
presents us with methods, which go beyond simple scal- 
ing level arguments, to tackle periodic sysiems with dis- 
order. These techniques were used to showtZl that despite 
pinning, a classical three dimensional periodic system is 
defect free and retains quasi long range translational or- 
der (Bragg glass). In this paper, we use these methods 
to obtain a quantitative theory of the pinned WC. We 
show that an elastic description is an excellent approxi- 
mation. We compute both the static and the a.c. trans- 
port properties and investigate in detail the magnetic 
field dependence of the various relevant quantities, such 
as the pinning frequency and threshold field. We show 
that the behavior of a pinned WC has marked differences 
from the one that was commonly accepted, based on the 
analogies with pinned charge density waves. Our theory 
permits a reasonably consistent description of the exist- 
ing experimental results, and suggests additional ways to 
check for a pinned WC through the transport properties. 
A summary of the method and some of the results on 
the a.c. transport were presented in a shorter form in 
Ref. 

The plan of the paper is as follows. In Sec. II, we define 
the elastic Hamiltonian used to study this system, and 
discuss the relevant length-scales. The Caussian varia- 
tional method (GVM) that we use to solve the problem, 
and the variational equations are described in Sec. 111. 
In Sec. [V, we obtain the solution of these equations. 



two sections and jump directly to Sec.^. In this sec- 
tion, we discuss the static properties of the pinned WC, 
i.e. the compressibility and the translational order. The 
dynamics is examined in Sec. 0. General features of 
the conductivity are described, together with the mag- 
netic field and density dependence of the a.c. transport. 
Our theory is comp ared with existing experiments on the 
pinned WC in Sec. VII, while we discuss other theories 
in Sec. 



VII] . Conclusions can be found in Sec. 



IX and 



many of the technical details have been relegated to the 
Appendices. 



II. MODELING THE SYSTEM 

A. Elastic model 

Let us define the model we use to study the insulat- 
ing crystal phase. Firstly, we see that retaining the full 
fermionic nature of the wavefunction, would lead to a 
very complicated problem in the presence of disorder and 
coulomb interactions. Fortunately, in the crystal phase 
the particles become localized, and thus distinguishable. 
So a "classical" description for the crystal is a good start- 
ing point. The quantum aspects of the original fermionic 
problem will be hidden in the parameters characteriz- 
ing this crystal. In the crystalline phase, the electrons 
occupy the sites of a triangular lattice with a lattice con- 
stant a which is related to the density of electrons by 
n ^ (7ra^)~^ . As shown in Figure ^ the electrons at site 
i are displaced from their mean equilibrium positions la- 
beled by the two dimensional vector i?^, by u(Ri, t). The 
"effective" size of the particles in this crystal are deter- 
mined by the spatial extension of the electronic wave- 
function, as shown in Figure 0. In the presence of a 
strong magnetic field, this is simply the cyclotron radius 
Ic = \f\J(eW). At finite temperatures, thermal fluctua- 
tions add to this purely quantum extension of the wave- 
function. Hence the "size" of the particle gets renormal- 
ized. For high temperatures, it tends to the Lindemann 
length /-r, defined by the average thermal displacement 



A rough estimate of the size of the particle 



These two sections contain most of the technical parts 
of the paper, and we concentrate on the physical results 
in the remaining sections. The reader interested only in 
the physical consequences of our theory can skip these 



is thus I'x- In the following, we will simply de- 

note this (temperature dependent) length by Ic- Note 
that this kind of classical approach is certainly valid for 
strong magnetic fields where, Ic ^ a. In this case, the 
overlap of the various electronic wavefunctions is negli- 
gible and the Fermi statistics can be safely ignored. It 
is rather difficult to give a precise estimate of the field 
and density for which this approach breaks down. A use- 
ful hint is nonetheless given by the Lindemann criterion 
for the melting of a classical crystal ~ C'£a'^, where 
the Lindemann number Cl ~ 0.1. This shows that even 
displacements close to melting remain relatively small. 
It is thus reasonable to expect that our approach would 
remain valid even relatively close to melting. Our ap- 
proach also applies to Wigner crystal a,t B — 0, provided 
the correct extension of the wavefunction is used in place 
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FIG. 1: The three length characterizing the Wigner crystal. 
The size Ic of the "particles" in the crystal (at low tempera- 
ture it is essentially given by the extension of the wavefunction 
around the equilibrium position and is the cyclotron radius, 
at large temperatures it is controlled by the thermal fluctu- 
ations and is the Lindemann length), a the lattice spacing 
is controlled by the density of particles, and the disorder is 
correlated over a length rf. The inset shows the triangular 
structure of the Wigner crystal. Particles are labeled by an 
equilibrium position Ri and a displacement it. 



of Ic- 

The long wavelength properties of the WC can now be 
described by an elastic hamiltonian describing the dis- 
placements u. We set ?i = 1 in the rest of the paper. 
The vibration modes of this crystal can be quantized in 
the usual manner leading to the imaginary time elastic 
action in the presence of disorder, 



+ iprnLO„U}c^al3\U(3{-q, ~UJn) 

+ [ (fr [ dTV{r)p{r,T) 



(1) 



where here and below we denote the integration over the 
Brillouin zone (BZ) 



BZ 



(27r)2 



(2) 



and Uq, with a ~ x,y are the two components of the 
vector u, eap is the antisymmetric tensor txy — —^yx = 1 
and the Matsubara frequencies Un = with the inverse 



temperature (3 = 1/T. p,„ — -^^Pc — ^ffi are the mass 
and charge densities, c and d are the shear and bulk 
moduli. These moduli can be obtained from an expansion 
of the coulomb correlation energy of the WC in terms of 
the displacements. Analytical expressions for the elastic 
moduli were obtained for the classical crystal in Rcf. |2|. 
Later, efforts were made to estimate quantupa. corrections 
to these moduli using trial wave functionsEj. Here, we 
work with the classical values obtained in Ref. where 
the elastic moduli are independent of the magnetic field, 
which should be a reasonable approximation for strong 
fields. The bulk modulus d and the shear modulus c are 
given by 



PI 
e 



Pla 



(3) 
(4) 



where a ~ 0.03 and e is the dielectric constant of the sub- 
strate. Since the bulk modulus describes compressional 
modes, it is drastically affected by the coulomb repul- 
sion. The off-diagonal term in the action (0) comes from 
the Lorentz force, and is proportional to the cyclotron 
frequency Wc = PcB/ pm, where B is the magnetic field 
applied perpendicular to the plane of the crystal. Finally, 
the last term describes the coupling to disorder, modeled 
here by a random potential V . The density of particles 



p{r) ^(5(r - i?i - Ui) 



(5) 



where (5 is a ^-like function of range Ic (see Figure |l|) and 
Ui = u{Ri). Since the disorder can vary at a lengthscale 
rf a priori shorter or comparable to the lattice spacing a, 
the continuum limit Ui — *■ u(r), valid in the elastic limit 
\ui — p4fepM I ^ a should be taken with care in the disorder 
termO'cJ. This can be done using the decomposition of 
the density in terms of its Fourier components 



(6) 



where po is the average density and K are the reciprocal 
lattice vectorspof the perfect crystal. The finite range of 
S is recoveredtZI by restricting the sum over K to mo- 
mentum of order iCmax t^/Ic We assume a gaussian 
distribution for the disorder V{r), which leads to the fol- 
lowing disorder averages 



V{r)V{r') = Arf{r~r') 



(7) 



is a delta-like function of range r/ which is the char- 
acteristic correlation length of the disorder potential (see 
Figure |l|). The gaussian limit is valid when there are 
many weak pins. We will come back to this point in a 
more quantitative manner in the next section. 
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B. Physical lengthscales 

The physical properties of the various phases of the two 
dimensional electron crystal are completely governed by 
the interplay of the the various microscopic length scales 
in the system: the lattice spacing a, the particle size 
Ic and the disorder correlation length rf. These micro- 
scopic parameters define collective lengthscales which are 
of fundament al irnpg rtance in the physics of such disor- 
dered systemsllIpEj'E3. These length scales, Rc and Ra, 
are the distances over which the relative displacements 
are of the order of the size of the particle and the lattice 
spacing respectively. They are defined by the disorder 
averaged displacement correlation functions 



(8) 



In order to distinguish between these two lengths it is 
mandatory to keep all the harmonics in the decomposi- 
tion of the density (^). Keeping only one harmonic, as is 
done for CDW, amounts to considering that the effective 
"size" of the particle is a and hence Rc ~ Ra- For the 
Wigner crystal, as we will see, it is crucial to carefully 
distinguish between these two lengthscales. 

The length Ra can easily be obtained by a static scaling 
argument comparing the cost in shear elastic energy and 
the gain due to disorder. For the WC, where n ~ (7ra^)~^ 
we get 



Ra 



(9) 



In order to be in the weak pinning regime, and for an 
elastic description to be valid, it is necessary to have 
-Ra > a. 

The interpretation of Rc, the so called Larkin- 
Ovchinikov lengthllil is more subtle. It corresponds to 
lengths where the particles can truly feel the random po- 
tential. In particular, it is the length above which the 



Larkin mode0, which approximates the random poten- 
tial by a random force becomes invalid, and pinning and 
metastability appear. If Rc 3> a, the system is collec- 
tively pinned, whereas one would have single particle pin- 
ning in the opposite case. We emphasize that the elastic 
approximation is valid regardless of the behavior of Rc, 
provided one still has Ra 3> a. 

Since u{Ra) ~ a, Ra has in the past been incorrectly 
interpreted as the lengthscale at which topological defects 
such as dislocations appear and the translational order of 
the crystal is lost. The naive picture is the one of a crys- 
tal broken into crystallites of size Ra- It has recently 
been shown that this view is grossly incorrect. In par- 
ticular, even in two dimension topological defects appear 
at a lengthscale Rd much greater than Ra, and the iCTjys.- 
tal preserves a very good degree of positional orderOtS. 
In three dimensions, the situation is even more favor- 
able since the system preserves a quasi translational or- 
der (power-law divergent Bragg peaks) and topological 
defects are not generated by weak disorder {Rd — oo). 
The physical implications of these results for the present 



problem will be discussed in section VI D 



III. GVM 



Method 



We now study the disordered WC described by (|^). 
Due to the nonlinear coupling of the disorder to the 
displacement field u in (0), this problem is extremely 
difficult, t|0 .solve. Here, we treat it using a variational 
methocOEJ. Many of the technical details and subtleties 
of the method can be found in the literature (see e.g. 
Rcf. ^^|J2C| for a review). We focus directly on the WC 
and we present here only the main steps. The first step 
involves averaging (|l|) over disorder by introducing repli- 
cas. This averaging results in an effective action which 
involves interactions between the n replicas, given by 



<(-9,-W„) 



^ / 1^ drdr' V Ak cos[K ■ {u''ir,T) - u\r,T'))] 
2 J Jo Jo „ I. ,^ 



(10) 



a,b,K 



where summations over a and /3 are implicit. The replica 
indices, a, b run from 1 to n. K denotes the reciprocal 
lattice vectors. The size of the particles Ic and the finite 
correlation length of the disorder restrict the sum over K 
to values of K smaller than i^max ~ tt/ max(rj, ^c)- Ak 
is taken as a constant Ak = A for K < iiTmax and zero 



otherwise. A more precise formulation is given in Ap- 
pendix The physical disorder averages are recovered 
in the limit n — *■ 0. (^) is quite general and can be used 
to descrjiie a host of physical systems, both quantum and 
classicalcJ. 
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The disorder couples the various rephcas via the co- 
sine term. Though this cosine term is local in space, it 
is completely non-local in time. This non-locality leads 
to a rich dynamical behavior as we shall show in this pa- 
per. In addition to the cosine term, another quadratic 
term is generated by the q = part of the disorder. This 
can be absorbed in a shift of the displacement field u. It 
affects the static correlation functions but not the trans- 
port properties and thus we neglect it henceforth. 

We now search for a variational solution to (J10| ) by 
using the best quadratic action approximating (|10D. We 
use the trial action 



(11) 

are vari- 
ational parameters. The variational free energy is now 
given by 



where the whole Green's function ^(9,^ 



{S-S, 



0)So 



(12) 



The variational parameters are then determined by the 
saddle point equations 



dGiy{q,cu,,) 



(13) 



The explicit form of these equations are given in Ap- 
pendix PI 



B. Saddle point equations 

The saddle point equations ( p^ have to be solved in 
the limit of the number of replicas n — *■ 0. There are 
two kinds of saddle point solutions in general: one with 
replica symmetry and the other with replica symmetry 
breaking. For the case of interest here, d = 2, one can 
show that the correct way tOr-take the limit n ^ is 
to break the replica symmetrytS. We first introduce the 
displacement correlation function 



Bf^{x,T) = (K(x,r)~u^(0,0)]2) = ^ ^(GrM«,'^«) + G^''/3(9,^n)-2cos(ga; + c.„r)GS^^(q,c^„)) (14) 

J q n 



r 



We then parametrize the replica Green's function as fol- 
lows: 



) = r^iq,iJr,)Sat+<j:,^M (15) 



a/3/ 



where a are the variational parameters and f°'^{q,uJn) 
is the elastic matrix defined in (Q). Since the disor- 
der induced interaction between replicas is local in space 
and non-local in time, the parameters a depend only on 
the frequency. We now define the connected part of the 
Green's functions as 



b 



(16) 



The saddle point equations, the fap and the various con- 
nected Green's functions are all given in the appendix 

The matrix structure of the saddle point equations sim- 
plify if we use the basis of longitudinal and transverse dis- 
placements. This offers a physically transparent picture 
in that the transverse modes are the shear modes of the 
solid and the longitudinal ones represent the compres- 
sion modes. In this basis, the Cartesian displacements 
are given by 



Uaiq) = u^{q)qa + {q)€ai3qp 



(17) 



where q — q/\q\ is the unit vector along q. For the pure 
system, in the absence of a magnetic field, the longitudi- 
nal and transverse modes are the eigenmodes of the sys- 
tem. The longitudinal mode is sensitive to the Coulomb 



repulsion whereas the transverse mode is not. The non- 
disordered part of the Hamiltonian can be rewritten as 



S[u] 



(PmW, 

q n 

U„ ,„ (Pmt^n + Cq )U_g 
T 



cq^ + dq)u^g_^ 



q,UJri 



-U)ri 



118) 



The magnetic field mixes these two modes. Note that 
the longitudinal mode yields the correct dispersion for 
the plasmon mode 



q 



1/2 



Though the disorder term 
in ( |10|) has a complicated form in terms of ul and mt, the 
saddle point equations themselves can be easily written in 
terms of these longitudinal and transverse modes. Since, 
the appropriate solution is the replica symmetry broken 
one, all quantities which are off-diagonal in the replica 
indices are parametrized in terms of a continuous variable 
< u < 1- There is a break-point Uc such that for 
u < Uc the solutions manifest explicit replica symmetry 
breaking and one recovers effectively replica symmetric 
solutions for u > Uc- The saddle point equations, their 
parametrisation in terms of u, the structure of the various 
replica Green's functions are all discussed elaborately in 
Appendix ^ 

In this section, we present only the final version of 
these equations. These equations contain a parameter E 
which can be seen as a disorder induced mass gap and 
the function liiujn) which describes essentially the dissi- 
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pation due to disorder. In terms of these new variables, the various connected Green's functions are now given by 
I 



GcL {q, jw„) 

GcLTiq,iuJn) 



cq^ + dq + Pm'^'i + I{iuJn) + S(l - 5nfi) 



[W - 






Cq^ + prr 


5nfi)){cq'^ + Pmi^n ^ 
UjI + I{iuJn) + S(l 


f I{iuJn) - 




- 5n.0)) - 


- 0^ up" uP'\ 


[(cg2 - 


\-dq + Pmi^l 






<5ri,o))(cg2 + Pm^l - 
PmLOnt^c 


f I{iuJn) - 




- ^n,o)) - 


^ pL^I^I] 


[(cg2 - 


\-dq + Prn(^l 




+ - 


5nfi)){cq'^ + Pmi^n ~ 


f I{iuJn) - 




- Sn.o)) - 


1- Pm'^n'^c] 



(19) 



where I{iuj„) satisfies (in the semi-classical limit) 



I{iuJn) 



1 



cq^ + dq 
2{cq^ + 



dq 



[cq^ + pm^^'i +dq + I{iujn) + S)(c(j2 + p^W^ _|_ /(i^j^J + S) + pla^'i^l 

I 



(20) 



and the equation determining S is 
if 

The term in the exponential, B(uc), is given by 
B{u,) = ■ 



(21) 



{GcL{q,^n) + GcT{q,^n)) 
1 1 



cg2 + (ig + E cq^ + E 
where the breakpoint is 



(22) 



(23) 



(determined in Appendix Since ( pi[ ) depends on / 
through B(uc), (jll), (|^) and (|2|) form a closed set of 
self consistent equations. 



IV. SOLUTION OF SADDLE POINT 
EQUATIONS 

We now solve the self consistent equations derived in 
the previous section. Before, we proceed with the explicit 
determination of the solution, let us examine some of the 
general features of the solution. A cursory glance at the 
propagator (p^, shows that the mass term E defines a 
lengthscale through E = d~^. As explained in Ref. |l^. 



this is the length scale which separates the Larkin regime 
from the random manifold regime. These two regimes are 
characterised by a power law growth of displacements 
with differing exponents. Moreover, at zero tempera- 
ture, in the absence of thermal fluctuations the relative 
displacement correlation function, Bt{1) oc max[/ 



Comparing with (| 
equivalcntly. 



we immediately see that I — Rc or 
E = cR^^ (24) 



where, Rc is the Larkin-Ovchinikov length. Contrary to 
older theories of CDW, E is related to Rc and not the 
positional length Ra- We emphasize that E = cR~^ only 
in the collective pinning regime where Rc 3> a. The 
equivalence / = Rc fails when the system enters the single 
particle pinning regime where Rc < a. 



1. Results for S and Rc 

To solve for E, it is useful to remark that the struc- 
ture of the equation ( ^l| ) , is very similar to the analogous 
equation for a classical system with point-like disorder. 
This allows us to directly apply the formalism developed 
in RcfJlJ to the present problem. Re-arranging the terms 
in (B19), we can see that B{uc) is the extent of local- 
ization of the particle. Therefore, to leading order, we 
replace B in (^ll) by which gives 



c{27:^)-R-\a/^o)' 



(25) 
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where 



^0 = max[r/,y 



The Larkin length is given by 

Rc = Rai — )^ 

a 



(26) 



(27) 



We see from (^5|), that E exhibits very different magnetic 
field dependences in the two physically distinct cases of 
r/ >< Ic- We see that for r/ Zc, S is independent of 
the magnetic field and for r/ <C Ic, S — bB^, where 
reinstating h, b = c{2tt'^)- ^ R'"^ {a"^ e / Tif (at T = 0). 
This has extremely interesting ramifications for experi- 
mentally measurable quantities as we discuss in the fol- 
lowing sections. 

Let us illustrate the importance of retaining all the har- 
monics in the case where disorder varies at length scales 
smaller than a. If we retain only the lowest harmonic in 
( pl| ) , we would have rf = a since disorder now varies only 
on scales comparable to a. In the relevant limit, Ic < a, 
we see that the exponential in (pl| ) can be replaced by 
unity. In other words, this corresponds to substituting 
^ = r^: = a in ( p5| ) leading to the result that Rc = Ra- 
This precludes the possibility of any dependence of S on 
the magnetic field because Ra is a static scale indepen- 
dent of B. Thus we see that when disorder varies at 
scales smaller than a i.e., r/ < a, it is necessary to retain 
all the harmonics in (|2|) to obtain the correct result for 



2. Behavior of I{uj) 

Though we have solved the equation determining E in 
the full quantum sense, ( B17 ) proves difficult to solve. 
We thus examine the semi-cl^ical limit given by (20). 
As was shown in simpler casesH, this limit captures cor- 
rectly the main features of the solution. (pO|) is analyt- 



ically solvable for small and large frequencies and can 
be solved numerically for intermediate frequencies. Two 
frequency scales control the behavior of /(u;„): (i) the 
pinning frequency ujp = ^/pmUJc, to be discussed in de- 
tail in section VI; (ii) the cyclotron frequency lUc- For 
small ujn (wn < j^r) ' S&^^ 



I{iuJn) 



'2p„,E 



7rp?„^g(cE)^ 
2d 



(28) 



As expected I(iujn) when w„ ^ 0. In the limit of 
small LUc, the second term inside the square root can be 
nelegected and we obtain the usual result for a system 
in a zero magnetic field. On the other hand, the second 
term dominates the low frequency behavior of I, if 



(29) 
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FIG. 2: A generic numerical solution for the real and 
imaginary l2{u>) parts of I{u>) (in units of S) as a function of 
Lo (in units of y/'E/pm)- Results for two different values of u)c 
are shown. The corresponding value of ujp = T,/pm<-Uc is also 
given. 



This defines a frequency scale, ^ . We work in the 
limit of large magnetic fields, Pm^^c ^ ^- Since the two 
elastic moduli are related c? oc -, fi^ cx where m* is 
the effective mass of the system. This clearly shows that 
is just the plasma frequency of the two dimensional 
WC. 



For intermediate frequencies 
solution becomes 



the 



V n2 , ,2, ,2 

6 °^ rfEl 



(30) 



The analytical solutions given above, were obtained in 
the limit of Pm^c ^ ^- To obtain the solution numeri- 
cally, we first continue ( pO| ) to real frequencies w, which 
results in a set of two coupled integral equations for the 
real and imaginary parts Ii and I2 of /(w). These in- 
tegral equations can be easily solved numerically. The 
numerical solutions indicate that the linear solution (|2^ ) 
obtained for frequencies uj <C Y^/ Pm^c, infact survives 
beyond this regime. The numerical approach, however, 
is not constrained to this regime and can be used to ob- 
tain / for all ranges of E and lOc- Typical plots of Ii and 
I2 are shown on Figure |^. We see from Figj^, that I2 
is indeed linear in lo and Ii is quadratic for small uj, in 
agreement of the analytical continuation of (p8|) . 
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V. PHYSICAL PROPERTIES - STATICS 
A. Compressibility 

One of the few thermodynamic properties that can be 
measured in the QHE and MOSFET systems is the static 
compressibihty k of the Wigner crystal. The formahsm 
used in this paper allows us to calculate the static com- 
pressibility of the disordered crystal. The calculation of 
K requires a knowledge of the density-density correlation 
function. 

K = lim lim I{{q,uj) (31) 

g— ►O >0 

where the density operator is given by daUa = QaUa- 
Clearly, only the longitudinal components of the displace- 
ment field contribute to the density. Consequently, the 
density correlation is entirely determined by the correla- 
tion of the longitudinal displacements and is given by 

n(g,c^) = plq'^ {uL{q,oj)uL{-q, -uj)) (32) 

First let us consider the pure case. We immediately see 
that since the coulomb term dq dominates over the shear 
term for small q, the compressibility k = in the pure 
system. Consequently, the compressibility remains zero 
even in the disordered phase. We reiterate that the com- 
pressibility calculated here is related to the second deriva- 
tive of the energy with respect to density for fixed total 



number of particles. This also assumes that the neutral- 
ising background remains unchanged as the volume of the 
electron gas is changed. Before comparing these results 
with experiments, it should be ascertained whether the 
compressibility measuijed in the experiment is the same 
quantity defined aboveo. 



B. Translational order 

It is interesting to study the effect of disorder on the 
translational long range order present in the clean sys- 
tem at zero temperature. Contrary to thcr^aive view 
based on extrapolation of Larkint2l and FLEO methods 
that the crystal loses its long range order due to pin- 
ning, it was shown that in c? = 3 the disordered system, 
although pinned, retains a defect-free quasi-long range 
translational order. In c? = 3, this Bragg-slass phase has 
thus power-law divergent Bragg peakalU'Ej. For d = 2 
classicaLsystems it is known that topological defects do 
appeared but at a length Rd which is much larger than Ra 
and thus the elastic description holds in a large regime of 
length scales (up to Rd) resulting in a quasi Bragg glass 
state. The properties of the d = 2 quasi Bragg glass for 
the classical problem have been studied-iii-detail and the 
length scale Rd has been estimatedllZrE§Eauj. 

We now derive the result given by the variational 
method for the positional correlation functions of the 
pinned WC in d = 2. 



2 

/3uc 



qdqdO 



duj 



(1 — cos(gxcos0)) 



GcL,t{GcL,T + S) 



(1 — cos^ 9){1 — cos{qx cos 9)) 



GcT,l{GcT,L + S) 



(33) 



Since the equal time asymptotic behavior is governed pri- 
marily by a; = 0, it is sufficient to retain this mode alone 
in (p3|). We find that the leading contribution to the 
correlation function grows logarithmically i.e., 

Bl..t{^) =^ T^-^ log ( 1 (34) 

Substituting (|23| ) in the above and using the relation S ~ 
cR~^ at zero temperature, ( p4[ ) simplifies to 

Only the leading contribution to the displacements grow 
in an isotropic manner. There are non-isotropic correc- 
tions to this result but these can be neglected in the 
asymptotic limit. The translational correlation is given 



by 

Ot = (exp^-^o-^^^) exp-*^° (36) 
Within the Gaussian approximation used in this paper 

KoaKo/S ~ 

Ot = cxp[ — -Baf3[x)\ 

thus the pinned WC also supports quasi long range order. 
Since Rc is determined by (p7|), we see from (|4|) that the 
growth of displacements will have a magnetic field depen- 
dence depending on whether is lesser or greater than 
Ic- The B dependence in a static correlation function like 
the one studied here is a direct consequence of the fact 
that in the quantum system the statics and the dynamics 
are coupled. 
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Although the GVM correctly captures several features 
of the growth of displacements in d = 2, it should be 
taken with a grain of salt concerning the exact asymp- 
totic behavior. It is well known for the classical problem 
that the variational method slightly underestimates the 
fluctuations at the lower critical dimension d = 2 and 
that the variational result log(r) should in fact be re- 
placed by log^(r). The latter result can be dcEivedJaj^ 
more accurate renormalization group approachOllj'E3E^. 
It is fair to expect that the result will be similar for the 
quantum problem considered here, since the tu = mode 
considered in ( ^ ) mimics a d — 2 classical problem. De- 
spite this, even a log^ growth of the correlation functions 
would lead to an extremely weak destruction of the posi- 
tional order, in the quantum case analogous to the quasi 
Bragg glass. The other important question is the role of 
dislocations. Again the precise answer is only known for 
the classical problem, for which it has been shown that in 
d = 2 dislocation are generated but at a lengthscale Rjj 
much largeiO than Ra- In particular the result ^atpzero 
temperature, relevant for the ground state, read£j'Lj: 



given by 



(38) 



where co is proportional to the strength of disorder. For 
the pinned WC the effect of additional quantum fluctu- 
ations and Coulomb interactions remains to be treated 
in a precise way. As can be seen e.g. by analogy to 
a 3d classical system with correlated disorder (where it 
is known that fluctuations play a minor role), one can 
surmise that (38) should be rather stable to small quan- 
tum fluctuations. Consequently, the pinned WC should 
still exhibit an extremely good positional order for weak 
disorder. 



VI. PHYSICAL PROPERTIES - DYNAMICS 
A. General features of the conductivity 

The formalism and the results obtained in the previ- 
ous sections can now be used to calculate the magneto- 
conductivities of the pinned crystal. The conductivities 
are given by 



(Ta0{(^) = iplujGap{q = 0, W + id) 



(39) 



a, (3 = x,y and the G are the displacement Green's func- 
tion. The real frequencies co are analytic continuations of 
the Matsubara frequencies a;„. Using (B3), the explicit 
form of the conductivities are 

2 iUj[~PmOJ^ + S + 



^(S-p„C^2+J(^))2_p2^^2^2 
2 Oj[pm<^(Mc] 



'yx 



(E - p,„tj2 + /(w))2 - pl^uluj 



^40) 



The conductivity tensor is therefore, completely deter- 
mined by S and I{uj). As usual the Hall resistivity is 



R 



H 



(41) 



'xy 



For the pure system, E = / = and (|40|) gives back 
the standard features of electrons in a magnetic fleld. 
The Drude peak is shifted by the magnetic field from 
u! — to the cyclotron frequency ujc- The d.c. Hall 
conductivity has a finite value crxy{uj = 0) = Pc/B, and 
the Hall resistivity has the classical value Rh = Rl Pc- 

In the presence of disorder (S ^ 0, /(w) 7^ 0), as can 
be seen from (^ ) , the (5- function peak at lOc will broaden 
and slightly shift in frequency. (|o|) also shows that the 
dc conductivities Oxxi^jJ = 0) = '^yyiy = 0) = cFxyiy — 
0) = 0, both in the presence and absence of a magnetic 
field. This is a consequence of the crystal being pinned by 
disorder. More importantly, a new peak appears in the 
diagonal conductivity at new scale LOp called the pinning 
frequency. 

A simple estimate of the pinning frequency is given 
by the poles in the conductivity when setting I{uj) = 0. 
This gives a. lo^ — {^/ujf+^Sjp^ — uJc)/'^- In the hmit 
of large magnetic fields, which is relevant for the case of 
quantum Hall systems, the above expression simplifies to 



(42) 



Of course a;° is 



The cyclotron peak is shifted to uJc + 
only an approximation to the true pinning frequency LUp 
that should be determined directly by the maximum of 
the conductivity (^. The main difference between ujp 
and LOp stems from the presence of /(w) in the denomina- 
tor. Typical values of the pinning frequencies are shown 
in Fig. I and Fig. |. 

Quite remarkably, the structure (40) of the conductiv- 



ity tensor shows that the Hall resistivity Rh is a constant 
independent of co 



Rh 



B 

Pc 



(43) 



as in the pure system. This is a remarkable result since 
since it shows that the disorder has no influence on the 
Hall resistivity of a pinned WC. This property comes 
from the fact that the disorder is local, which in the 
variational calculation makes all off diagonal self ener- 
gies zero, as explained in Appendix Since this feature 
is linked to the locality of disorder it is expected to be 
valid beyond the variational approximation. 

We now analyze the frequency dependence of the con- 
ductivities. For small values of u), substituting ( p8|) in 
{M gives 



Reaxxiio) 
Reaxyico) 



Pc\ "^Pni^ 



2d 



S2 



2 \2 
PcPmi^ci^) 



(44) 
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° FIG. 5: r/ > Ic'. Keaxx (in units of pc/\/pmS) as a function 

of u! and a;c(in units of -^E/pm). Note the non Lorentzian 
FIG. 3: rf > I,: The pinning frequencies ujp and as a gj^^pg ^^e peaks. In this regime, the pinning frequency 
function oi luc- All frequencies are in units of ^S/p™,. Note decreases with increasing magnetic field, 
that Up and cop become quite different for low magnetic fields. 





FIG. 4: ff < Ic- The pinning frequencies lUp and ijjp (in 
units of pc/pmb) as a function o f the applied field B (in units 
of Pc/pmb). b is defined in Sec. [V 1, 



FIG. 6: rf > Ic'. lmaxx,^eaxy and Imaa;j;(in units of 
Pc/Vpm'S) as a function of uj and uic (in units of y^T./pm) 
for Lo^ = 10. 



In the region 



< w < — we find using (30) 



and the imaginary or dissipative part of the conductivity 
grows as 



2 ^ 



(45) 



(46) 



The full behavior of the conductivity can be obtained 
from ( ^ ) by numerically solving the equation for I{u)). 
Typical plots of diagonal and transverse conductivities 
are shown in Fig. ^ and Fig. ^. 

In the limit of very strong magnetic fields, since uj'^ ^ 

<^ ujc, provided VI J ^ (see discussion in 
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FIG. 7: Keaxx and Imaxx (in units of p^/y' pm'S) as a func- 
tion of (jj{m units of y^T,/pm) for the pinned crystal in zero 
magnetic field. 



Sec. [V 2) the evolution of the width Aujp of the pin- 



ning peak is completely governed by the low frequency 
behavior of / obtained in (Eq). In this limit, 



cx - 

ujp \n 



1/4 



(47) 



Again from (|2^), we see that the behavior of the peak 
widths is again dictated by and Ic- A similar result 
was obtained in Ref. ^ with an undetermined exponent 
s (to be discussed in Sec. VIII). 



The case of zero external magnetic field can also be 
obtained from ( ^0|) by setting = 0. It leads to different 
energy scales and physical behavior. The clean crystal 
has infinite dc conductivity and zero ac conductivity. In 
the presence of disorder, the dc conductivity is zero as 
expected and a pinning peak develops at the frequency 
ujp = ^Hj pm where S is determined by (|2|) with = rf. 
We show a typical plot of the conductivity in Fig. 0. 

An important feature of our results is that the pinning 
peaks are rather sharp for finite magnetic fields, as op- 
posed to the rather broad pinning peak obtained in the 
limit of zero magnetic field. This is expected because the 
cyclotron peak in the former case has appreciable spec- 
tral weight. Comparing Fig. |^ and Fig. |, we see that 
the peaks in the finite field case are at least a factor of 
10-15 times narrower. 

Moreover, the older theories predict a Lorentzian 
broadening of the peaks. A glance at the conductivities 
plotted in Figs||,^ and ^ shows that the pinning peaks 
are non-Lorentzian and asymmetric around LOp. This is 
an important prediction of our theory. This asymme- 
try in the peak is indeed seen in microwave. pQpductivity 
measurements on certain Ga-As samplesEj'oEj. Com- 



parison of our results w ith e xperiments will be discussed 
in further detail in Sec. VII. 



A generic property of pinned systems is that a finite 
nonzero threshold electric field Et is required be- 
fore the crystal de-pins and d.c. current flows through 
the system. Since Et is intrinsically the signal of non- 
linear response, it cannot be directly evaluated using 
Kubo-liJ«|J3rmulae. As for the dynamics of classical 
crystals'^ •iri, its precise calculation would require a non 
equilibrium formulation such as the Keldysh technique. 
This is clearly a formidable task, but fortunately for the 
classical cases, approximate derivations allow one to re- 
late Et to equilibrium quantities. A dynamical deriva- 
tion has shown that such formulas indeed yield the cor- 
rect threshold field Et- Here, we use the same approxi- 
mate approach to derive Et for the quantum system. In 
the collective pinning regime, a, the threshold field 

Et is determined by i?JiJ 



Et - cR-'^o 



(48) 



where = niax[rf , Ic]. Clearly, Et defined by ( p^ ) is 
directly linked to ojp by Et — pmi^c(^p£,a- The validity of 
this result remains 1;o be explicitly proven in the quantum 
case. 



B. Magnetic Field dependence 

As already pointed out, it is crucial to have theoretical 
predictions for the various features and attributes of the 
conductivity as a function of physical parameters like the 
magnetic field and the density. This would then facilitate 
direct comparisons between experiment and theory and 
consequently further our understanding of the system. 

It is clear from the previous section that E and hence 
the length scale Rc determine most of the dynamical 
properties of the pinned crystal. Different behaviors are 
possible depending on the relative magnitude of the cor- 
relation length of the disorder ry and the particle size Ic 
(see Fig. |l|). This leads to a physics very different from, 
that predicted by the naive weak pinning scenarioEZIEj 
where, it is implicitly assumed that the scale Ra is the 
only lengthscale in the system. Two main behaviors can 
be obtained. 



> Ic 



In this case, the particle locally sees a smooth poten- 
tial, and thus behaves effectively like a point particle. 
Hence the length Ic is not directly relevant. This leads to 
a Rc which is independent of the magnetic field, as was 
shown in Sec. IV. ^From (Eq), one obtains for the naive 



pinning frequency 



-1 cR„ 



PcB Tf 



(49) 
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Consequently, the naive pinning peak moves towards the 
origin as and gets narrower with increasing field. 

However, as mentioned in Sec. VIA, the true pinning 



frequency is given by the location of the peak in the op- 
tical conductivity, and can be in principle different from 
LUp. A plot of the conductivity in this regime for different 
values of the magnetic field is shown in Fig. |^. From it we 
extract the true pinning frequency ujp. As seen in Fig. 
both frequencies decrease as B^^. We see that the differ- 
ence ujp — Ljp is more pronounced for small values of B or 
equivalently, w^/S, and goes to zero for increasing values 
of the field B. Similarly, the non-Lorentzian nature of the 
peaks discussed earlier is enhanced for small and inter- 
mediate magnetic fields (cf. Fig. ||) . The regime rf > Ic 
includes the case of the conventional CDW in a magnetic 
field for which rf ^ a and thus Ra — Rc- Therefore, 
CDW have the same magnetic field dependence for the 
pinning frequency as the Wigner crystal in this regime. 
The density dependen ces ar e, however, quite different, as 
will be shown in Sees. VIC and VI] . 



2. Tf < Ic 

In this regime, since the size of the particle is larger 
than the correlation length of the disorder (see Fig. |^), 
the particle locally averages over the disorder. The ef- 
fective disorder seen by the particle thus strongly de- 
pends on its size l^.. Since the particle size is strongly 
magnetic field dependent, this regime exhibits anoma- 
lous magnetic field dependence in all physical quantities. 
This physics is captured by the equations for S . (|25| ) 
indicates that S is strongly B dependent. This explicit 
magnetic field dependence results in a naive pinning fre- 
quency (we have reinstated fi) 



cR- 



Pc 



B^ ,ea^„ 



(50) 



which increases quadratically with increasing field B. 
Therefore, unlike the previous case ( ^9|) where the pin- 
ning frequency and peak width were decreasing functions 
of i?, here the pinning peak moves up and broadens with 
increasing field. 

it is important to extract 



As was done in Sec. VI B 1 



the true pinning frequency and peak width from the full 
conductivity. Rea^x has been plotted for various values 
of B ranging intermediate fields to strong fields, in Fig. ||. 
In this regime, all curves are plotted for varying w^/^- 
However, since E cx B^, ^cf^ ^ B~^, which implies 
that an increase in w^/S corresponds to a decrease in the 
magnetic field B. A plot of the true pinning frequency is 
shown in Fig. ^. One sees that in this case the shift in 
pinning frequency, compared to the naive expression ( |50| ) 
is quite important. As before, this difference is more pro- 
nounced for small and intermediate values of w^/E which 
contrary to the previous case, now corresponds to large 
values of the field B. In a similar vein, as seen in Fig.(@), 
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FIG. 8: rf < 1^: Kecrxxiin units of bpm/pc) as a function ofuj 
(in units of p'^/p'^b; b has been defined in tlie text). Note that 
increasing <x;^/E corresponds to decreasing magnetic fields. 



the peaks become more and more non-Lorentzian for in- 
creasing B i.e., decreasing uj^/T,. Our result for the field 
dependence of LUp, shown in Fig. ^, can be compared with 
the experimental results. This point will be discussed in 
further details in Sec. VIl. 

As long as rf < Ic, uJp increases with B. This increase 
in the pinning frequency is not indefinite. Two effects 
can limit this regime. For sufficiently large B, since Ic 



decr eases w ith B one will cross over to the regime of 
Sec. VI B 1 for which Ic < rf. Since from (|5|) and (|27 



Rc decreases with increasing B, this regime is also limited 
by Rc ~ a, below which the system is no longer in the 
collective pinning regime. In that case, the correspon- 
dence (E4|) between Rc and S no longer holds and the 
result (pSF) for E is invalid. One should thus re-evaluate 
E in this single particle pinning regime. Solving (Kj) in 

3 

this regime, we find E cx i? 2 and 



uj° cx B^ 



(51) 



Which of the two crossovers {Ic ~ rf or Rc ~ a) occurs 
first depends on the strength of the disorder in the system 
considered. 



3. Summary of magnetic field dependence 

In general, one can thus expect three regimes with dra- 
matically different behaviors of experimentally measur- 
able quantities. A summary of these regimes and their 
features is shown in Table Contrary to the standard 
CDW case, for a Wigner crystal the pinning frequency 
can increase with the magnetic field provided one is in 
the regime rf < Ic- In all the regimes, the height of the 
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Regime 


E 




Et 




rf > Ic 


B" 


B-' 


B° 


B° 


rj < Ic 




B^ 


si 


B-' 


Rc < a 


Bi 


Bi 


B 


B-i 


n (r/ >< Ic) 


1 


3 


1 




n ^ 


n 2 


n 2 





TABLE I: Magnetic field and density dependences of various 
dynamical quantities. 



pinning peak decreases approximately as B^^ with in- 
creasing field. As expected, the cyclotron peak at LOc + tOp 
is not qualitatively affected by the relative magnitudes of 
rf and Ic and always moves upwards with increasing B. 

The difference between and ujp that appears in the 
WC regime, urges a re-examination of the relation be- 
tween the pinning frequency and the threshold field. In- 
deed if the formula (^) holds, the threshold field should 
be related to and not ujp. This would imply that 
in the WC regime, the pinning frequency and threshold 
field would not be linearly related. It would thus be in- 
teresting to have extended plots of Et vs the pinning 



frequency Wp, such as the ones of Ref. g. Let us reiterate 
that although the pinning frequency ajp can be system- 
atically calculated, the expression ( [48| ) for the threshold 
field remains unverified in the quantum case. Assuming 
(^), the magnetic field dependence of the the threshold 
field Et and the dielectric constant defined by 



ImUxx/^ 



(52) 



are summarized in Table. 



C. Density dependence 

Another way to probe the Wigner crystalline state is 
to fix the magnetic field and see how the pinning peak 
evolves as the density of electrons is varied slowly. The 
first thing to realize is that since r/ and Ic are both in- 
dependent of the density, the density dependence of S in 
( p5| ) and hence lo^ are the same irrespective of the rel- 
ative mag nitudes of and Ic- The pinning frequency 
using (p5|) is given by 



eoA 



muc e^^l 



(53) 



( p3| ) shows that the pinning frequency increases as den- 
sity is decreased. The width of the peak increases and 
the height of the pinning peak decreases as -n? with de- 
creasing n. This is analogous to the behavior seen for in- 
creasing B when rj < Ic- Accordingly, other dynamical 
quantities like Et, Sxx etc have the density dependences 
indicated in Table. || 



D. Limitations and dc conductivity 

Before embarking on a comparison of our results with 
experiments, let us discuss the range of validity of our 
theory. 

As already mentioned in Sec. II A, the quantum na- 
ture of the original particle can be safely hidden in the 
parameters of the crystal hamiltonian, as long as the spa- 
tial extent of the particles remain small compared to the 
lattice spacing. One expects this to be valid even quite 
close to the melting transition of the WC. Of course, our 
approach does not allow for a very precise description of 
the melting properties, beyond the simple position of the 
transition, since this simplified description is incapable 
of taking into account the fermionic statistics and the in- 
compressible quantum fluid that is the melted "crystal" . 

An a priori limitation of our study, even in the in- 
sulating phase far from melting, is to know whether the 
elastic model is a satisfactory starting point to obtain the 
transport properties or whether defects (vacancies, inter- 
stitials, dislocations etc) should be taken into account 
from the start. Indeed if dislocations appear at a typi- 
cal lengthscale Rd they will strongly modify the results 
given by the elastic description for frequencies smaller 
than a characteristic frequency w^i oc R]j^. 

Older theories of disordered elastic systems suggest 
that due to the presence of disorder, defects are favored 
at the length Ra for which displacements would become 
of order of the lattice spacing (see section II B ) . If, as 
in FL description, the pinning peak would be associated 
with the length Ra, we see that dislocations would change 
most of the interesting optical conductivity and in par- 
ticular the pinning peak. Such a situation is summarized 
in Figure §(a). If this was correct this would seriously 
impair any description of the transport properties based 
on an elastic hamiltonian. 

In fact, both of these commonly accepted points are 
incorrect, which makes the description based on an elas- 
tic model much more useful than initially anticipated. 
Firstly, as we already discussed in Sec. VI A| the pinning 
peak is not associated with the length Ra but with the 
a priori much smaller lengthscale Rc (note that larger 
Icngthscales correspond to smaller frequencies). Sec- 
ondly, as discussed in Section VB defects in disordered 



elastic systems are much less important than initially 
believed and appear at a lengthscale Rjj much larger 
than Ra {Rd = oo in d = 3 for weak disorder). Since 
Rd > Ra > Rc, the frequency lut) is thus much smaller 
than the characteristic frequencies of the pinning peak 
LUp. This implies that dislocations will not drastically af- 
fect the pinning peak and thus that the whole pinning 
peak structure is given quantitatively by the elastic 
model. This is summarized on Figure |^(b). Our study 
is thus perfectly adapted to describe optical conductivity 
experiments. 

For d.c. transport or static properties the role played 
by topological defects is less obvious. In particular, it 
is reasonable to assume that they can have an impact 
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VII. EXPERIMENTS 





FIG. 9: (a) If dislocations occurred at scale Ra and the pin- 
ning frequency was controlled by Ra, as was naively believed, 
the elastic theory is incapable of giving any reliable informa- 
tion on the pinning peak. It would be necessary to include 
dislocations from the start, (b) As was shown in the text, 
dislocations occur in fact at Rd S> Ra and the pinning peak 
depends on Rc <^ Ra- Thus the pinning peak is given quanti- 
tatively by a purely elastic theory. It is necessary to take into 
account topological defects such as dislocations only at much 
lower frequencies, and in particular if one wants to obtain 
reliable results for the d.c. transport. 



on some of the static properties. Their relevance for the 
thermodynamics is hard to estimate but very Hkely to 
be reasonably small since the contribution of the elastic 
part of the hamiltonian is usually finite and thus dom- 
inant. A notable exception might be the compressibil- 
ity. However, dislocations are mostly expected to play 
a role in the d.c. transport properties. Indeed a pinned 
system, in the absence of topological defects has a zero 
d.c. li near r esponse to an external force, as was shown 
in Sec. VIA. The calculation of the nonlinear response 



has only been done for classical systems for which it has 
been shown that to move, the system has to de-pin col- 
lectively larger and larger bundles, as the external force 
F becomes small, leading to a very small velocity of the 
form V ~ e^Pi'^/P)'^ _ For a quantum elastic systep. it is 
yet unclear how this creep law should be modifiedES, and 
what would be the d.c. transport of the purely elastic sys- 
tem. Nonetheless, the existence of topological defects al- 
lows for a linear, even if exponentially small, response to 
exist, and this should dominate the d.c. linear response. 
The temperature dependence of the d.c. resistivity is also 
clearly linked to the nature and mobility of the defects 
in the crystal. Clearly, obtaining the dynamics of a dis- 
ordered system in the presence of topological defects is a 
formidable task. This again emphasizes the importance 
of optical conductivity experiments to analyse systems 
expected to be pinned elastic structures, since they per- 
mit a comparison with controlled theories. 



Several experiments measured the optical properties of 
a two dimensional electron gas subjected to a strong mag- 
netic field. One important question is of course whether 
these data, obtained in the regime of small filling frac- 
tions where crystallization is expected, can be satisfac- 
torily explained within a pinned WC theory, or whether 
more exotic descriptions are necesaaii^. 

A first batch of experimentstj'ESES obtained the pin- 
ning frequency and the general frequency profile of the 
pinning peak for various values of the applied magnetic 
field. Important features of these experiments were a 
pinning peak whose frequency increased with increasing 
magnetic field. This qualitative behavior is in contradic- 
tion with the n|aive expectation based on the Fukuyama- 
Lec argunicntsEZI, that the pinning frequency decreases 
with increasing B. This was interpreted as casting doubts 
on whether such peaks could be interpreted within a 
pinned WC phase. As we pointed out, this conclusion 
is incorrect, and the experiments are in fact consistent, 
within our theory, with a pinned WC. Indeed, as we 
pointed out the FL-type theories incorrectly assume that 
Ra is the lengthscale controlling the pinning peak, while 
the relevant lengthscale is Rc- 



As was discussed in Sec. IV 2, the detailed dynami 



cal properties depend on the disorder and on whether 
^ riy^Ti/pm or not. The typical experimental pa- 
rameters for the hole-like samples are: B 13T, Uc = 
7.6 X lO^^Hz, m* — 0.3me, where nie is the bare elec- 
tron mass, e — 13eo (eo is the dielectric constant of the 
vacuum) and densities of the order of n = 5 x lO^^cm"^. 
For these values, using the classical values of the elastic 
moduli, the plasma frequency = 4 x Since 
the measured pinning frequencies are of the order of 
u)p ~ 7 X lO^Hz, we see that the condition ( p9[ ) is easily 
satisfied, implying that the experimental systems lie in 
the strong field regime discussed in Sec. VIA. Here, the 



pinrpag frequen cy can in fact increase with the magnetic 
fielcO (cf. Sec. VIA). Thus the observed behavior is in 
fact compatible with the one of a pinned WC. Our the- 
ory also predicts a non lorentzian shape for the peaks, as 
showa. in Fig. ||, a feature that is clearly exhibited by the 
dataB. 

However several features seen in these experiments are 
still extremely puzzling, and in fact incompatible, to the 
best of our knowledge, with all existing theories. Some 
anomalous features are the saturation of the pinning fre- 
quency and the increase of spectral weight in the pinning 
peak when the magnetic field in increased. Within a WC 
approach, if one in addition assumes field independent 
elastic constants, one always expects a reduction of the 
spectral weight contained in the pinning peak due to the 
continuous transfer of spectral weight from the pinning 
peak at low frequencies to the cyclotron peak. On the 
other hand, it is crucial to note that the interpretation 
of the experimental data is complicated by the following 
facts: (i) the experiment is performed at a quite high 
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temperature T = bQmK which is comparable with the 
pinning frequency of / ~ IGHz. Therefore, thermal 
contributions to the broadening of the peaks might have 
to be taken in. (ii) The filling factors considered are quite 
high. It is thus plausible that most of the data lies out- 
side the range of validity of a putative WC. It would thus 
be extremely useful, to have more detailed experiments 
performed in a lower range of temperatures and at lower 
filling fractions, (iii) Since the experiment is performed 
at relatively large filling factors, the very hypothesis of 
field independent elastic constants, which is valid at large 
fields, is questionable. One would have to know the true 
dependence of the elastic constants on the magnetic field, 
to make a reliable comparison even within the framework 
of a pinned WC theory. 

Another series of experimentslll focused on the den- 
sity dependence of the pinning peak. This set of experi- 
ments were done at lower temperatures T — 25mK and 
lower filling fractions, so this data is much more likely to 
be in a regime of a WC. Since these experiments are at 
fixed magnetic field, they partially avoid the above men- 
tioned problem with the elastic constants. They should 
thus provide for a more clear cut interpretation. A de- 
crease of the pinning frequency with increasing density 
was observed, fitting reasonably well with either a n^^ 
or a n~2 behavior, depending on the sample. A com- 
parison of these results was made with the FL theory for 
CDW. Within this theory and for a CDW, the pinning 
frequency varies as n~2. A blind application of this re- 
sult to the WC led Ref. |l^ to conclude that there was 
good agreement between the theory and the experimen- 
tal findings. But as shown in Appendix the correct 
extension of the FL theory to the WC case, leads in fact 
to an increase of the pining frequency to density as . 
This is due to the fact that Refs. |2^,^ deals with charge 
density waves where the equivalent of the displacement 
field u is the phase which is dimensionless. Reinstating 
the correct factor a which is the unit of displacement 
gives an extra factor of the density n^. The fact that a 
naive extrapolation of the FL theory is not possible for 
pinned WC is again clearly demonstrated by the experi- 
ment. The correct behavior ( p3| ) is n~2. So our theory 
is perfectly compatible with the experimental data. In 
addition, other features such as the non-lorentzian shape 
of the peaks, are also clearly visible in the data. In this 
experiment, the peak height decreases and peak width 
increases as the pinning frequency increases. This is in 
qualitative agreement with o ur predictions for the pinned 
WC, as shown in Sec. VIC. Compare, for example the 



qualitative lineshape of Fig. g with the data in Fig. 2 of 
Ref. HJ. Thus the agreement of these experiments with 
our theory strongly supports the interpretation of the 
observed phase as a pinned WC. 

Furthermore, our results suggests other possible exper- 
iments to check for further experimental signatures of the 
WC. The middle frequency range (i.e. Wp ^ w ^ Uc) 
could be compared with the theoretical formula (46). 
Working at frequencies above the pinning peak has two 



advantages: (i) the system is more free from thermal 
fluctuations, since frequency becomes the dominant en- 
ergy scale; (ii) as explained in Sec. VI D the effects of 
the topological defects are totally unimportant, so one 
expects our variational results to be a very good descrip- 
tion of this regime. Another useful probe of the existence 
of the crystal, is given by difference in behavior of the 
conductivity as a function of the both the density and 
the magnetic field. In the liquid phase, all quantities are 
functions of the filling fraction f = n/B, but do not de- 
pend separately on these two quantities. The striking 
difference between the n and B dependence seen in ex- 
periment is thus another strong indication of the presence 
of the crystal. It would thus be interesting to systemat- 
ically investigate whether a scaling in n/B is obeyed or 
not. Such a scaling could then signal the melting of the 
crystal. 



VIII. COMPARISON WITH OTHER THEORIES 

Let us now compare our work with other theoretical 
work on the conductivity of the Wigner crystal. 

As already mentioned in the previous sections, our 
work and the pioneering work of FL for the pinning of 
CDW share the same starting point. However, it differs 
from the FL-type analysis in the following essential ways: 
(i) the GVM permits a complete and consistent calcula- 
tion of the dynamical conductivity. This allows for a 
precise calculation of the pinning frequency and the dis- 
order induced dissipation (see Sec. VI B 3| ) ; (ii) the GVM 
bypasses the problem of the unknown fitting parameter 
that occurs in the Self Consistent Born ApproximationnZl, 
and thus provides an unambiguous determination of the 
conductivity. The GVM allows for a proper description 
of the inherent glassy aspects of the system, and in par- 
ticular, of the metastable states; (iii) more importantly, 
to treat correctly the case of the pinned WC, it is nec- 
essary, contrarily to the CDW case, to retain variations 
of disorder at lengthscales smaller than the particle size. 
This point was overlooked in more recent applications of 
the FL theory to the case of the pinned Wdl^O. As 
has been extensively discussed in the preceding sections, 
this leads to radically different dependence of the pinning 
frequency on the magnetic field and density. 

The question of the linewidth of the pinning peak in 
conductivity was also addressed in two series of works. 
Refs. p9pO| map the the crystal onto an effective bosonic 
model. Within a mean field approach, they obtained ex- 
ponentially narrow widths of the pinning peak. Since the 
computed width is much narrower than observed, extra- 
neous modes of broadening need to be invoked to explain 
the experimental results. In Ref. |3^, an expression for 
o'xx was obtained in the limit of pinning frequencies much 



smaller than the cyclotron frequency. From the conduc- 
tivity expression a formula, linking the half width of the 
peaks with the low frequency part of the conductivity 
was obtained giving Aujp/ujp ~ [ojploJcY , where s is the 
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exponent governing the low frequency part of the con- 
ductivity. Using a quite different method, the authors of 
Ref. ^ unknowingly rederived formulas for the conduc- 
tivity (their formulas (62) and (68) and the associated 
energy scales) that are identical the ones obtained by us 
via the GVM in Ref. |l|. Not surprisingly, Ref. || find 
the same energy scales. The only difference is the dissi- 
pation function /(w) which is replaced in Ref. ^ by an 
arbitrary dissipation function (denoted / and taken to 
be purely imaginary) scaling for low frequencies as uP'^ . 
In our case the GVM completely fixes the function /(w) 
(with both a real and imaginary part) which scales as w^, 
thus fixing s = 1/2. In Ref. |3^, though the exponent s 
could not be computed, arguments were given that the 
value should be s = 3/2, leading to much more narrow 
peaks than those stemming from the GVM. Although, 
it is in principle possible for the GVM to miss a part of 
the excitations (see e.g. the discussion in Ref. p^|20 ) we 
argue that for the pinned WC the exponent s = 1/2 is 
correct for the following reasons: (i) while it is conceiv- 
able for the GVM to miss part of the excitations (such 
as solitons) this oversight should lead to an underesti- 
mate of the low frequency part of the conductivity, while 
an exponent s = 3/2 would imply that the GVM over- 
estimates the conductivity; (ii) in d = 1, the exponent 
s = 1/2 given by the GVM, is exactly the one needed to 
reproduce the correct behavioij-Qf the conductivity which 
can also be determined exactlylla. In addition the param- 
eter s = 3/2 leads to peak widths that are much smaller 
than the ones that are observed experimentally. Perhaps 
the difference between our results and the approach fol- 
lowed in Ref. |3| to estimate s, can be traced to the fact 
that Ref. |3^ uses essentially a classical approximation for 
the WC. The /2(aj) oc w on the other hand depends on 
?i so is a quantum effect and will vanish if all quantum 
effects are thrown away. 

In addition to these approaches based on an elastic 
description of the system, Ref. ^incorporated a different 
approach based on the microscopic model of an electron 
gas confined to its lowest landau level. Though they find 
that the pinning frequency increases with the field, they 
cannot account for the sharp increase in peak amplitudes 
nor do they obtain the correct dispersion of the collective 
modes. The authors suggest the interesting possibihty 
that electronic single particle excitations might in some 
indirect way be responsible for the experimentally seen 
behavior. Further work is required before unambiguous 
theoretical conclusions can be drawn. 



IX. CONCLUSIONS 



pinning peak, can be reliably extracted from an elastic 
theory, without having to take into account topological 
defects of the pinned WC. We have studied the evolu- 
tion of various dynamical quantities such as the pinning 
frequency when the magnetic field and density are var- 
ied. We have shown that this problem is controlled by 
quite different lengthscales as compared to the standard 
pinned CDW problem. More specifically, depending on 
the strength of the magnetic field, the pinning frequency 
can exhibit quite different behaviors as a function of mag- 
netic field. Our theory presents a consistent explanation 
of most of the observed expjgcknental features, particu- 
larly the recent experimental3'Ej'E3 and thus provides a 
strong support for the interpretation of the high field 
phase observed in the 2DEG as a pinned Wigner feffH 
tal. Note that the deviations in the Hall coefficientcSH 
observed in the pinned phase, are c onsis tent with the pre- 
dictions of a transverse critical forcecScj that should exist 
for a pinned periodic structure, again in agreement with 
the interpretation of this phase as a pinned Wigner crys- 
tal. Conductivity and threshold measurements on the 
zero-field WC phase seen in Ga-As samples, as a func- 
tion of B and n should provide an ideal testing ground 
for our theory. Also a better understanding of the disor- 
der present in these samples is required before a realistic 
comparison with theoretical models can be made. 

Our work prompts several questions and issues that 
still need to be analyzed. An important issue is of course 
the influence of the anomalous (fractional quantum hall) 
liquid phase. As we have argued this should have little 
influence deep in the crystal phase, but certainly changes 
the melting or the properties close to melting in a way 
that still remains to be understood. This might be a 
reason for the strange behavior observed at intermedi- 
ate magnetic fields. The way the disorder or the elastic 
moduli could be affected by the magnetic field in this 
region also demands attention. But the most important 
question is the effect of topological defects at finite tem- 
perature. As we have shown, they are not important for 
the a.c. transport. But most likely they will be central 
in understanding the d.c. transport. It is not known 
whether the dc resistivity has any anomalous temper- 
ature dependence like those seen in experiments where 
Pxx{T) exp{To/TY where v is an exponent which is 
close to 0.5 in Si MOSFets and is found to vary with 
density in the experiments on the zero field GaAS sam- 
ples. Other interesting questions concern the scaling with 
temperature and electric field seen in the data. Can such 
a scaling be accommodated within the pinned WC ap- 
proach ? We leave these directions for future work. 



We have presented in this paper a quantitative theory 
of the pinned Wigner crystal. We computed static and 
a.c. transport properties. We have obtained the complete 
frequency dependence of the imaginary and real parts of 
the conductivity tensor. We have shown that the a.c. 
transport properties and in particular, the features of the 
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APPENDIX A: EFFECTIVE DISORDER 

Here, we derive the disorder term present in ( [lo| ) . Av- 
eraging (|^) over disorder leads to an inter-replica inter- 
action of the form 



Sint = -y drdr' J dridr2^rf{r - r')6{r - ri) 

S{r' - r2)pi^=Q{ri,Ti)pi^=o{r2,T2) (Al) 
(the time integrations are implicit) where 

P = ^ <^(r - i?^ - u,) (A2) 

i 

is the density for point like particles. Using (with an 
unbounded summation over K) one easily gets 



5'int = -plof ^ J drdr' J dridr2Arf{r - r')S{ 



r - ri) 



(A3) 



Given the nature of the various functions in (A3) and 
that ri — r2 ^ a (typically ri — r2 < max(r/, Ic)) , in the 
elastic approximation, one can replace u(ri, ti)— M(r2, T2) 
by u(r,Ti) — u(r,T2) wh ere r = (ri + r2)/2. Then one 
easily obtains (|l^) from ( |A3| ) with 



Ak = A{K)5{K)5{~K) 



(A4) 



where A{K) and 5{K) are the respective Fourier trans- 



forms of and 5. 



APPENDIX B: SOLUTION OF THE 
VARIATIONAL EQUATIONS 



Here, we present the derivation of the saddle point 
equations (|20| ) and ( pl| ) and the detailed structure of the 
variational solutions. Using the parametrization of the 
replica Green's function given in ([l5|), we minimize the 
free energy to obtain the following saddle point equa- 
tions. The off diagonal elements of a yield 

alf ^-^Jdr cos{u;^t)V^p {B'^^') (B1) 

where the displacement correlation function B is given 

by (ii 



V{B 



E 

K 



A 



(B2) 



and the prime denotes derivatives of V with respect to 
the argument B^^^ . The diagonal components yield an- 
other set of equations involving the connected part of the 



Greens functions defined by G^^p = 



Gcu 

Gc22 
Gcl2 
Gc21 



1/3 



{Y + Fn)iX + F22) -{Z + Fi2)(-Z + ^^21) 
{Y + Fi^){X + F22) -{Z + Fi2){-Z + 



Z + Fu 



(Y + Fii){X + F22) -{Z + Fi2){-Z + F21) 

-Z + F21 

(y + Fii){X + F22) ~{Z + Fi2){-Z + F21) 

where 1 and 2 denote the coordinates x and y. 



X = (p^cul + cq' + d^) 
Y = (p^iol + cq^ + d^) 

Z = {praLOc^n+ d'^^^) 



(B4) 



and 



Fo^p = ^ j dT{l-cos{uj.nT))Y,^KK^K, 



K 



(B5) 



The above equations form a closed and self-consistent set. 
These have to solved in the limit where the number of 
replicas n — > 0. Typically, there are two generic solutions: 
one is the replica symmetric solution and the other is the 
replica symmetry broken solution. The solution appro- 
priate to the problem studied depends on the number of 
space dimensions, the regime of temperature studied and 
the nature of the disorder. For the present problem of 
a weakly pinned WC, it has been showntS that the rele- 
vant solution is the one with replica symmetry breaking 
(RSB). Accordingly, we parametrize the replica Green's 
functions in terms of hierarchical matrices such that 



Gap{q,i^n) 
Gap{q,i^n,u) 



(B6) 
(B7) 



where the continuous variable < u < 1. The dis- 
placement correlations B given by ( |l4[ ) are similarly 
parametrised, resulting in two functions Bap{u) and i?a/3- 
An important simplification occurs when one consid- 
ers quantum problems with quenched disorder or equiv- 
alently classical systems with correlated disorder: since 
in each realization of the random potential, the disor- 
der is time independent, all quantities off-diagpaal in the 
replica indices are always r— independentiianj. Using 
this property, (Bl) now simplifies to 



(B8) 



Thus all off-diagonal replica coupling, and in particular 
the replica symmetry breaking (RSB) is confined entirely 
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to the LUn = mode. For the present problem, an addi- 
tional simplification accrues from the fact that (B2) in- 
volves only the local (i.e. at r = 0) S and B. These local 
quantities are isotropic and thus Baf3{0,u) = 5a[3B{u) 
and BajsiOjT) = 6af}B{T). This isotropy results in the 
enormous simplification 



As a result, 



Fi2 — F21 — 



(B9) 



(BIO) 



Therefore, the numerator of Gci2 in ( |B3| ) is unaffected 
by disorder. 

We now define a new variable [a] — ua{u) — a{v)dv. 
Based on the generic form of the replica solutions, we 
introduce the break-point u^. such that cr(u) = cr{uc) for 
u > Uc- If we retain only the lowest harmonic in the 
cosine term, the solution has one step RSB with a{u) — 
for u < Uc and <j{u) = cr(uc) for u > Uc- In terms of the 
variable [a], this RSB ansatz takes the form 



[fT](w) =0 u < Uc 
[a]{u) = = 2(3ucV' {B) 



U> Uc 



Physically, the region u < Uc corresponds to the loga- 
rithmic regime, characterized by a logarithmic growth of 
displacements and m > Wc is the Larkin regime where 
[cr](w) = E for u > Uc- This is basically a regime where 
for length scales below the Larkin length Rc, the physics 
of the system can be approximated by that where inde- 
pendent random forces act on each electron. However, for 
a multi-cosine model there exists a third kind of behav- 
ior for intermediate length scales. This regime, called the 
random manifold regime, is characterized by a power law 
growth of displacements, but with an exponent v which 
is different from that in the Larkin regime. 

The above discussion can be summarized in technical 
terms as follows: the solution [a] = is valid upto a 
break-point Ua where the random manifold regime starts 
manifesting itself. This is characterized by a [a] oc 
where /i is the random manifold exponent. This power 
law growth of [a] is valid upto a second break-point Uc 
beyond which the Larkin kind of behavior takes over. For 
single cosine models like CDW, Ua = Uc and the random 
manifold regime shrinks to zero. The absolute value of 
S and hence _Rc, depend on the the number of Fourier 
components of the disorder that is retained in the prob- 
lem. Hence it is important to retain all the harmonics in 
the density expansion so as to get the correct estimate of 
S. Here, for the quantities of our interest, though it is 
necessary to retain the higher harmonics in the expres- 
sion for E, we do not require the explicit value of Ua for 
our purposes and and hence it is sufficient to obtain the 
value of Uc alone. 

We now switch to the transverse and longitudinal com- 
ponents described in the main body of the paper. These 



render the saddle point equations more tractable and 
yield the following connected Green's functions 



g: 



{cq^ + p^ul) + F 



{cq^ + dq + p^Lul + F) 



(B12) 



Gel = [cq^ + dq + PraUjl) + F + 
1 



2 2 2 



(C(?2 + p„w2 + 

{cq^ + p,nLol + F){cq^ + dq + p,nOjl + F) 

Pm^nLOc 



Using ( Bll ), the disorder dependent part of ( B12 ) can 
always be written as 



F = I{iiUn) + E(l - 6n^o) 



(B13) 



Re-casting all the saddle point equations in terms of the 
longitudinal and transverse coordinates, we find that in 
the n — > limit .these reduce to equations determining 
Uc, E and /(iw„)cj. E is determined by 



1 = -2V"{Biuc)) f 



1 



1 



_(cg2-^S)2 {cq^ + dq + 



(Bll) where the off diagonal element B{uc) is given by: 



B{uc) 



+ 



^(Gcl((?, W„) -I- GcT{q,^n)) 
1 1 



cq^ + (ig + E cg2 + S 



(B15) 



and retaining only the fundamental harmonic ifo in 
we obtain 



-"■0 
87rc 



(B16) 



We draw attention to the fact that Uc is unaffected by 
the magnetic field and has the same value as that for a 
classical system with coulomb interactions and pointlike 
disorder. Lastly, 



/(ZC^„) = - / dT(\ - COs{uJnT))[V' {B{t)) - V'{B{Uc))] 

(B17) 



In ( B14 - B17 ), the local diagonal correlation 



B{t) = -((^(0,r)-u(0,0))2) 
^ \[Mr)+BT{r)] 



where 



Bt,l{t) = ^ [y^GcT.Li^ ~ cos(w„t)] 



(B18) 



(B19) 



In these expressions, GcL,t — (G^/t) ^- Since B{uc) 
depends on X^n^^o Gc(<^n) we see from ( p2| ) that E is a dy- 
namical quantity and hence can depend on the cyclotron 
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frequency Wc- To solve these equations, we work with 
a finite momentum cut off. Fortunately, for the case of 
the WC, the bulk modulus being much greater than the 
shear modulus, provides a natural cu toff A = -. 

Doing the momentum integral in (B14), the equation 
for S can be rewritten as: 



We then take the semi-classical limit Ti in (B17), 
This limit allows us to expand the exponentials in V. 
Using ( B16 ) and substituting the expressions for B{uc) 
and B 



K 



1 



K^B{u,) 



(B20) 



I{iujn) = -2V"{Q) 



I{iuj„ 



dq 



cq^ + S cq^ + dq + Y, {cq^ + Pm^^ + + I{iuJn) + S)(cg2 + p^w^ _|_ I(iujn) + S) 



r 



(B21) 



We have thus reduced the entire set of saddle point equa 



tions (|B1 



B3|) to two equations for S and I{ioJn) [ (|B20| ) 



and (B21)] which can then be solved self consistently to 



obtain the physical properties of interest. 



APPENDIX C: FUKUYAMA-LEE RESULTS FOR 
THE WC AND CDW 

In this appendix, we clarify the connection between the 
results presented here, and those of Ref. ^ as applied to 
the WC. The Fukuyama-Lee results have been misquoted 
in the literature. The pinning frequency of a CDW or WC 
obtained in Ref.E^ is given by 



The shear modulus c (x ni and Ra oc n^/^ for CDW 
and for the WC, Ra oc rt~2. This difference, stems from 
the fact that for CDW, the equivalent of the displacement 
field M is a phase variable, which is dimensionless. Within 
this approach, the density dependence is given Wpo ~ 
112 for the WC and ujpo ^ n~i for the CDW. This last 
result was erroneously used to compare the experimental 
data with a pinned WC in Ref. As our work has 
emphasized, the correct pinning frequency is determined 
by Rc and this leads to 



-3/2 



(C2) 



cR- 



(CI) 



for the WC. 
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